CONSERVATION OF GEOMETRIC STRUCTURES FOR 
NON-HOMOGENEOUS INVISCID INCOMPRESSIBLE FLUIDS 



Francesco Fanelli 



SISSA 

Functional Analysis and Applications Sector 
wo Bonomea 265 - 34136 Trieste, ITALY 



Universite Paris-Est 



Laboratoire d'Analyse et de Mathematiques Appliquees, UMR CNRS 8050 
61 Avenue du General De Gaulle - 94010 Greteil Gedex, FRANGE 



f rcLticesco . f cinelliOsissa. it 



f rcincesco . f EinelliSmath. cnrs . f r 



May 7, 2013 



Abstract 



We obtain a result about propagation of geometric properties for solutions of the non-homogeneous in- 
compressible Euler system in any dimension N > 2. In particular, we investigate conservation of striated 
and conormal regularity, which is a natural way of generalizing the 2-D structure of vortex patches. The 
results we get are only local in time, even in the dimension N = 2; however, we provide an explicit lower 
bound for the lifespan of the solution. In the case of physical dimension iV = 2 or 3, we investigate also 
propagation of Holder regularity in the interior of a bounded domain. 

1 Introduction 

In this paper we are interested in studying conservation of geometric properties for solutions of 
the density-dependent incompressible Euler system 

' dtp + u ■ V p = 



which describes the evolution of a non-homogeneous inviscid fluid with no body force acting on 
it, an assumption we will make throughout all this paper to simplify the presentation. Here, 
p{t, x) € R+ represents the density of the fluid, u{t, x) G its velocity field and Il{t, x) G M 
its pressure. The term VII can be also seen as the Lagrangian multiplier associated to the 
divergence-free constraint over the velocity. 

We will always suppose that the variable x belongs to the whole space . 

The problem of preserving geometric structures came out already in the homogeneous case, 
for which p = 1 and system ([1]) becomes 



in studying 2-dimensional vortex patches, that is to say the initial vorticity f^o is the characteristic 
function of a bounded domain Dq. As we will explain below, in the case of higher dimension > 3 
this notion is generalized by the properties of striated and conormal regularity. 



(1) 



< p {dtu + u ■ Vn) Vn = 
, div ti = , 



{E) 




divii = 



1 



The vorticity of the fluid is defined as the skew-symmetric matrix 



def:vort (2) 



n := Vn 



Vu 



and in the homogenous case it satisfies the equation 

dtn + u-vn + n-Vu + *Vn -n = 0. 

In dimension = 2 it can be identified with the scalar function u = div? — d2U^, while for 
= 3 with the vector-field uj = V x ti. Let us recall also that in the bidimensional case this 
quantity is transported by the velocity field: 



The notion of vortex patches was introduced in [26j and gained new interest after the survey 
paper |24) of Majda. In the case N = 2 Yudovich's theorem ensures the existence of a unique 
global solution of the homogeneous Euler system, which preserves the geometric structure: the 
vorticity remains the characteristic function of the evolution (by the flow associated to this solu- 
tion) of the domain Dq. Vortex patches in bounded domains of were also studied by Depauw 
(see |16)). while Dutrifoy in |17) focused on the case of domains in M'^. Moreover, in [6] Chemin 
proved that, if the initial domain has boundary ODq of class C^'*" for some e > 0, then this reg- 
ularity is preserved during the evolution for small times; in he also showed a global in time 
persistence issue. In Danchin considered instead the case in which initial data of the Euler 
system are vortex patches with singular boundary: he proved that if ODq is regular apart from a 
closed subset, then it remains regular for all times, apart from the closed subset transported by 
the flow associated to the solution. 

In the case > 3 one can't expect to have global results anymore, nor to preserve the initial 
vortex patch structure, because of the presence of the stretching term in the vorticity equation. 
Nevertheless, it's possible to introduce the definition of striated regularity, which generalizes in a 
quite natural way the previous one of vortex patch: it means that the vorticity is more regular 
along some fixed directions, given by a nondegenerate family of vector- fields (see definition 12.11 
below). This notion was introduced first by Bony in |3] in studying hyperbolic equations, and 
then adapted by Alinhac (see [Ij) and Chemin (see j^) for nonlinear partial differential equations. 

In [18j, Gamblin and Saint-Raymond proved that striated regularity is preserved during the 
evolution in any dimension N > 3, but, as already remarked, only locally in time (see also |25)). 
They also obtained global results if the initial data have other nice properties (e.g., if the initial 
velocity is axisymmetric) . 

As Euler system is, in a certain sense, a limit case of the Navier-Stokes system as the viscosity of 
the fluid goes to 0, it's interesting to study if there is also "convergence" of the geometric properties 
of the solutions. Recently Danchin proved results on striated regularity for the solutions of the 
Navier-Stokes system 



in |10) for the 2-dimensional case, in jl2) for the general one. Already in the former paper, he 
had to dismiss the vortex patch structure "stricto sensu" due to the presence of the viscous term, 
which comes out also in the vorticity equation and has a smoothing effect; however, he still got 
global in time results. Moreover, in both his works he had to handle with spaces of type -BpJ^ 
(with p g]2,-|-(X)[ and e €]2/p, 1[) due to technical reasons which come out with a viscous fluid. 
Let us immediately clarify that these problems have been recently solved by Hmidi in f21] (see 
also [2]). In the above mentioned works Danchin proved also a priori estimates for solutions of 
(NS,^) independent of the viscosity z^, therefore preservation of the geometric structures in passing 
from solutions of (NS,^) to solutions of (E) in the limit v ^ 0. 



dtUJ + u ■ Vw = . 




dtu + u-Vu - uAu + vn = 



2 



In this paper we come back to the inviscid case and we study the non-homogeneous incompress- 
ible Euler system ([T|). We will prove that preservation of geometric properties of initial data, such 
as striated and conormal regularity, still holds in this setting, as in the classical (homogeneous) 
one. 

Let us point out that considering geometric structures is not only a generalization of the 
classical problem of vortex patches. It can be also seen as an improvement to the well-posedness 
result for ([!]) in critical Besov spaces B^^. (see paper |14)). As a matter of fact, here we will 
make lower regularity hypothesis on the initial data, which will be compensated by the additional 
geometric assumptions. 

Note that in the 2-dimensional case the equation for the vorticity reads 

dtuj + u-Vuj + V A Vn = , 

so it's not better than in higher dimension due to the presence of the density term, which doesn't 
allow us to get conservation of Lebesgue norms. This is also the reason why it's not clear if 
Yudovich's theorem still holds true for non-homogeneous fluids: having ojq £ L*^ D combined 
with suitable hypothesis on po, doesn't give rise to a local solution. 

So, we will immediately focus on the general case N > 2. We will assume the initial velocity uq 
and the initial vorticity to be in some Lebesgue spaces, in order to assure the pressure term 
to belong to L^, a requirement we could not bypass. As a matter of fact, VII satisfies an elliptic 
equation with low regularity coefficient, 

- div (a vn) = div F , 

and it can be solved independently of a only in the energy space L^. Moreover, we will suppose Qq 
to have regularity properties of geometric type. Obviously, we will require some natural but quite 
general hypothesis also on the initial density pQ of the fluid: we suppose po to be bounded with 
its gradient and that it satisfies geometric assumptions analogous to those for Qq- Let us point 
out that proving the velocity field to be Lipschitz, which was the key part in the homogeneous 
case, works as in this setting: it relies on Biot-Savart law and it requires no further hypothesis on 
the density term. Let us also remark that no smallness conditions over the density are needed. 
Of course, we will get only local in time results. Moreover, we will see that geometric structures 
propagate also to the velocity field and to the pressure term. 

Our paper is organized in the following way. 

In the first part, we will recall basic facts about Euler system: some properties of the vorticity 
and how to associate a flow to the velocity field. In this section we will also give the definition of 
the geometric properties we are studying and we will state the main results we got about striated 
and conormal regularity. 

In section [31 we will explain the mathematical tools, from Fourier Analysis, we need to prove 
our claims: we will present the Littlewood-Paley decomposition and some techniques coming 
from paradifferential calculus. In particular, we will introduce the notion of par avector- field, as 
defined in [12\ : it will play a fundamental role in our analysis, because it is, in a certain sense, 
the principal part of the derivation operator along a fixed vector-field. Moreover, we will also 
quote some results about transport equations in Holder spaces and about elliptic equations in 
divergence form with low regularity coefficients. 

In section SI we will finally be able to tackle the proof of our result about striated regularity. 
First of all, we will state a priori estimates for suitable smooth solutions of the Euler system ([T]). 
Then from them we will get, in a quite classical way, the existence of a solution with the required 
properties: we will construct a sequence of regular solutions of system ([!]) with approximated 
data, and, using a compactness argument, we will show the convergence of this sequence to a "real" 
solution. Proving preservation of the geometric structure requires instead strong convergence in 



3 



rough spaces of type (for some a > 0). The uniqueness of the solution will follow from a 
stability result for our equations. In the following section, we will also give an estimate from 
below for the lifespan of the solution. 

Finally, we will spend a few words about conormal regularity: proving its propagation from 
the previous result is standard and can be done as in the homogenous setting. As a consequence, 
inspired by what done in Huang's paper [22J, in the physical case of space dimension = 2 or 3 
we can improve our result: we will also show that, if the initial data are Holder continuous in the 
interior of a suitably smooth bounded domain, the solution conserves this property during the 
time evolution, i.e. it is still Holder continuous in the interior of the domain transported by the 
flow. 



eq:BS-law (3) 



2 Basic definitions and main results 

Let u, VH) be a solution of the density-dependent incompressible Euler system ([T|) over [0, T] x 
and let us denote the vorticity of the fluid by Q,. As in the homogeneous case, it will play a 
fundamental role throughout all this paper, so let us spend a few words about it. 

From the definition (j2|), it is obvious that, for all q € [1, -|-cxd], if Vu € L^, then also € L'^. 
Conversely, if u is divergence-free, then for all 1 < i < A^ we have Att* = "^jLidjQij, and so, 
formally, 



N 



This is the Biot-Savart law, and it says that a divergence free vector-field u is completely deter- 
mined by its vorticity. From (|3]) we immediately get 



eq:BS_grad (4) 



N 



Now, as the symbol of the operator —di{—/S) ^ dj is a{^) = the classical Calderon- 

Zygmund theorem ensures tha10 for all g G ]1, +oo[, if Q £ L'^ then Vn e and 



est :CZ 



(5) 



\Vu\ 



Li 



<c- 



1 



Li 



In dimension N = 2 the vorticity equation is simpler than in the general case due to the 
absence of the stretching term. Nevertheless, as remarked above, the exterior product involving 
density and pressure terms makes it impossible to get conservation of Lebesgue norms, which was 
the fundamental issue to get global existence for the classical system (E). So, we immediately 
focus on the case N >2 whatever, in which the vorticity equation reads 



eq:vort (6) 



dt^ + u-vn + n-Vu + *Vii • V I - I A vn 



0. 



where, for two vector-fields v and w, we have set A if to be the skew-symmetric matrix with 
components 

{v A 'w)^j = v^iif — v^w^ . 

Finally, recall that we can associate a flow ip to the velocity field u of the fluid: it is defined 
by the relation 



ip{t,x) = ipt{x) ■= X 



u{t, ipr{x)) dr 



^This time the extreme values of q are not included. 
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for all {t, x) € [0, T] x and it is, for all fixed t G [0, T], a diffeomorpliism over M^, if Vu G L°°. 
Let us remark that the flow is still well-defined (in a generalized sense) even if u is only log- 
Lipschitz continuous, but it is no more a diffeomorpliism (see e.g. chapter 3 of [2] , or [8] , for more 
details). 

Let us now introduce the geometric properties we are handling throughout this paper. The 
first notion we are interested in is the striated regularity, that is to say initial data are more regular 
along some given directions. 

So, let us take a family X = (-'^A)i<A<m vector-fields with components and divergence 

of class for some fixed e G ]0, 1[ . We also suppose this family to be non-degenerate, i.e. 



I{X) := inf sup 



N-l , ^ 
A Xa(x) 



1 

JV-l 

> 0. 



Here A G Ajy_i means that A = (Ai, . . . , Aat-i), with each \ G {1, . . . , m} and Aj < Aj for i < j, 
while the symbol A Xa stands for the element of M such that 

V y G , {^K^ X^-Y = det (Xa, . . . Xx^_^,Y) . 

For each vector-field of this family we put 

IIXaIIc^ := IIXaIIc^ + lldivXAlIc^, 

while we will use the symbol |[| • [|[ in considering the supremum over all indices A G A™ = 
{1 . . . m}. 

d:stri Definition 2.1. Take a vector-field Y with components and divergence in and fix a 77 G 

[e, 1 -|- e]. A function / G L°° is said to be of class along Y , and we write / G Cy, if 
div (/y) G (R^). 

\i X = (^A)i<A<-m is a non-degenerate family of vector-fields as above, we define 

C\:= fl and ||/||^, := |7|X|||c. + ||| div (/X) 

l<A<m ^ ' 

r:div Remark 2.2. Our aim is to investigate Holder regularity of the derivation of / along a fixed 
vector-field (say) Y , i.e. the quantity 

2=1 

If / is only bounded, however, this expression has no meaning: this is why we decided to focus on 
div(/y), as done in the literature about this topic (see also |12) . section 1). Lemma 14.51 below 
will clarify the relation between these two quantities. 

Now, let us take a vector- field Xq and define its time evolution X{t): 
def:X| (7) X{t,x) = Xt{x) := dx,(cc)4^t {^rH^)) , 

that is X{t) is the vector-field Xq transported by the flow associated to u. From this definition, 
it immediately follows that [X{t) , dt + u ■'V] =0, i.e. X{t) satisfies the following system: 

f {dt + u-\/)X = dxu 



eq:X (8) 

We are now ready for stating our first result, on striated regularity. 



X\t=Q — Xo . 
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t:stri-N 



t : conorm-N 



Theorem 2.3. Fix e € ]0, 1[ and take a non- degenerate family of vector-fields Xq = (^o,A)i<;)^<m 
over R^, whose components and divergence are in . 

Let the initial velocity field uq € L'^, with p g]2,+oo], and its vorticity r^o € L°° D L^, with 
q e [2,+oo[ such that 1/p + 1/q > 1/2. 

Suppose the initial density pQ € W^'°° to be such that < p^, < < p* ■ 
Finally, let us assume that JIq (^i^d ^ PQ both belong to C^Xq- 

Then there exist a time T > and a unique solution (/), ti, VII) of system such that: 

• p G L°°{[0,T];W^'°°) nCb{[0,T] x R^), such that < p^ < p < p* at every time; 

• u £ C([0,r];LP)nL°°([0,r];C°'i), withdtu G C([0,r];L2) and vorticity Q G C([0, T]; L?); 

• Vn G C{[0,T]-L'^), withV^U G L°°([0,r];L°°). 

Moreover, the family of vector-fields transported by the flow still remains, at every time, non- 
degenerate and with components and divergence in , and striated regularity is preserved: for all 
t G [0, T], one has 



vp{t) , m G 



and 



uit) , vn(t) G c]+l^ 



uniformly on [0,T]. 



Another interesting notion, strictly related to the previous one, is that of conormal regularity. 
First of all, we have to recall a definition (see again section 1 of |12]). 

Definition 2.4. Let S C be a compact hypersurface of class C^'^. Let us denote by 7^ 
the set of all vector-fields X with components and divergence in C^, which are tangent to S, i.e. 
dxH |s = for all local equations H of S. 

Given a rj £ [e, 1 + e], we say that a function / G L°° belongs to the space Cv if 



\/x G r# 



div ifX) G C 



-1 



Similarly to what happens for striated regularity, also conormal structure propagates during 
the time evolution. 

Theorem 2.5. Fix e g]0, 1[ and take a C^'^ compact hypersurface Sq C M^. 

Let us suppose the initial velocity field uq G L^, with p G ]2, +oo], and its vorticity J^o G L°° fl L'^, 

with q G [2,-|-oo[ such that 1/p + 1/q > 1/2. 

Assume that the initial density pQ G Ty^'°° be such that < p^^ < po < p* . 
Finally, let Qq and V/Oq belong to . 

Then there exist a time T > and a unique solution (p, u, VH) of system ([1]), which verifies 
the same properties of theorem \2.'A 
Moreover, if we define 

S(t) := (So) , 

Ti{t) is, at every time t £ [0,T], a hypersurface of class C^'^ of M.^ , and conormal regularity is 
preserved: at every time t G [0, T], one has 



vp{t) , m G c| 



and 



uit) , vn(t) G 4+^ 



uniformly on [0,T]. 



3 Tools 

In this section we will introduce the main tools we used to prove our results; they are mostly 
based on Fourier analysis techniques. Unless otherwise specified, one can find the proof of all the 
results quoted here in [2j, chapter 2. 
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d:besov 



p:CZ 



3.1 Littlewood-Paley decomposition and Besov spaces 

Let us first define the so called "Littlewood-Paley decomposition", based on a non-homogeneous 
dyadic partition of unity with respect to the Fourier variable. So, fix a smooth radial function x 
supported in (say) the ball i?(0, |), equal to 1 in a neighborhood of .6(0, |) and such that r i— )■ 
x{t e) is nonincreasing over R_|_ for all unitary vector e G 

The dyadic blocks (Aj)jg2 are defined bjH 



if j < -2, 



-1 •- 



Moreover, set (^) = x {i/2)-x {i) ■ 
X{D) and := ip{2-W) if j > 0. 



We also introduce the following low frequency cut-off operator: 



E 

k<j-l 



Afc for j>0. 



The following classical properties will be used freely throughout the paper: 

• for any u (z S', the equality u = Aju holds true in <S'; 

• for all u and v in S', the sequence (Sj-iu Ajv)^^^ is spectrally supported in dyadic annuli. 
One can now define what a (non-homogeneous) Besov space Bp ^. is. 

Definition 3.1. Let u be a tempered distribution, s a real number, and 1 < p,r < +oo. We 
define the space Bp ^. as the set of distributions u (z S' such that 



\U\\B' 



< +00. 



From the above definition, it is easy to show that for all s € M, the Besov space S| 2 coincides 



s 

00,00 



is actually 



with the non-homogeneous Sobolev space H'^, while for all s G IR+\N, the space i?, 
the Holder space C^. 

If s G N, instead, we set := B^ ^, to distinguish it from the space of the differentiable 
functions with continuous partial derivatives up to the order s. Moreover, the strict inclusion 
C| ^ C| holds, where denotes the subset of functions bounded with all their derivatives 
up to the order s. 

If s < 0, we define the "negative Holder space" as the Besov space B^ ^. 
Finally, let us also point out that for any G N and p G [1, +00], we have the following chain 
of continuous embeddings: 

B'p^.^W'^'P^B^^^, 
where W^'^ denotes the set of functions with derivatives up to the order k in L^. 

Besov spaces have many nice properties which will be recalled throughout the paper whenever 
they are needed. For the time being, let us just mention that if the conditions 



s > 1 H 

P 



or 



1 -\- — and r 
P 

0,1 



1 



of bounded Lipschitz 



holds true, then Bp ^ is an algebra continuously embedded in the set C 
functions, and that the gradient operator maps Bp^^ in Bp~^. 

The following result will be also needed. 

Proposition 3.2. Let F : W be a smooth homogeneous function of degree m away from a 

neighborhood of the origin. 



Then for all {p,r) G [l,-|-oo] and all s G the operator F{D) maps Bp^ in B, 



s—m 
p,r 



Throughout we agree that f{D) stands for the pseudo-difTerential operator ut-^ T ^{f Tu) 
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1 :bern 



The following fundamental lemma describes, by the so-called Bernstein's inequalities, the way 
derivatives act on spectrally localized functions. 

Lemma 3.3. Let < r < R. A constant C exists so that, for any nonnegative integer k, any 
couple {p, q) in [1, +00]^ with 1 < p < q and any function u G L^, we have, for all A > 0, 

1 1 A 



suppuC B{0,XR) =^ llV'^nlli, < C^+1a''+^(^ '^II^^IIlp; 



LP 



c : embed 



As an immediate consequence of the first Bernstein's inequality, one gets the following em- 
bedding result. 

Corollary 3.4. The space Bp^^^_^ is continuously embedded in the space B^^ /o'" indices 
satisfying pi < p2 and 

S2 < si — N ( ) or S2 = si — N ( ) and ri < r2 . 

\Pi P2J \Pi P2J 

From now on we will focus on the particular case of Holder spaces. 
3.2 Paradifferential calculus 

Let us now introduce Bony's decomposition of the product of two tempered distrubutions u and 
v: we will define the paraproduct operator and recall a few nonlinear estimates in Holder spaces. 
Constructing the paraproduct operator relies on the observation that, formally, the product uv, 
may be decomposed into 



eq:bony (9) UV = TuV + TyU + R{u,v) , 

with 

TuV := ''^^ Sj-iu AjV and R{u,v) := AjuAj^v. 
j j \k-j\<l 

The above operator T is called "paraproduct", whereas R is called "remainder". 

The paraproduct and remainder operators have many nice continuity properties. The following 
ones will be of constant use in this paper. 



p:op 



c : op 



Proposition 3.5. For any s € M and t > 0, the paraproduct operator T maps L°° x into 
and C^* X into C^*, and the following estimates hold: 

WTuvWc < C ||^i||L°° llcs-i and \\Tuv\\(2s-t < C ||^i||c-i ||Vt;||cs-i . 

For any si and S2 in M such that si + S2 > 0, the remainder operator R maps x C'^ into 
QS-1+S2 continuously. 

Combining the above proposition with Bony's decomposition ([9|, we easily get the following 
"tame estimate": 

Corollary 3.6. Let u he a hounded function such that Vu € C^^^ for some s > 0. 

Then for any v £ we have uv £ and there exists a constant C, depending only on N 
and s, such that 

WuvWc < C^j|u||2,oo \\v\\cs + . 

In our computations we will often have to handle compositions between a paraproduct operator 
and a Fourier multiplier. The following lemma (see the proof e.g. in |10) ) provides us with 
estimates for the commutator operator. 
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est : comm 



p : comp 



l:coinni Lemma 3.7. Let m S M, i? > and f E C°°(IR^) be a homogeneous smooth function of degree 
m out of the ball B{0, R). 

Then, there exists a constant C , depending only on R, such that, for all s ^M. and all a <1, 
one has: 

C 

(10) \\[Tu, f{D)\v\\^s-m+a < ^ _ ^ \\Vu\\c^-i WvWc ■ 

Let us now quote another result (see [l^ for the proof of the former part, ilOj for the proof 
of the latter), pertaining to the composition of functions in Besov spaces, which will be of great 
importance in the sequel. 

Proposition 3.8. (i) Let I be an open interval o/M and F : / — )• M a smooth function. 

Then for all compact subset J C I and all s > 0, there exists a constant C such that, for all 
function u valued in J and with gradient in C^^^ , we have V{F o u) S C^^^ and 

\\V{Fou)\\cs^i < C|tVn||cs-i . 

(a) Let s > and m G N 6e such that m > s. Let u £ and tp G such that the Jacobian of 
is bounded. 

Thenuoip G . Moreover, i/s€]0,l[ the following estimate holds: 

\\u o ip\\^^ < C (1 + llVV'll^.oo) WuWc ■ 
Finally, let us introduce the notion of paravector-field. 



d:pvec-f Definition 3.9. Let X be a vector-field with coefficients in S'. We can formally define the 
paravector-field operator Tx in the following way: for all u S' , 



N 



Txu := ^ TxrdiU. 



1:T X 



i=l 

The following result (see |12]. section 2 for the proof) says that the paravector-field operator 
is, in a certain sense, the principal part of the derivation dx- 

Lemma 3.10. For all vector field X G and all u G C^, we have: 

• if t < 1 and s + t > 1, then 

(J 

\\dxu - Txu\\cs+t-i < ^-^ _ ^ _ ^-^ W^lh l|Vu||ct-i ; 

• if t < 0, s < 1 and s + t > 0, then 

t{s + tj'{s-l) 



Txu — div {uX)\\^s+t-i < — ;— ;^;-7 — \\X\\cs ll^llc* 



• if t < 1 and s + t > 0, then 

C ~ 

\\dxu - Txu\\(.s+t^i < ll-'^llc'' l|Vn||ct-i • 

Moreover, first and last inequalities are still true even in the case t = 1, provided that one replaces 
W'S/uWqo with llVnllioo, while the second is still true even if t = 0, with \\u\\loo instead of \\u\\qo. 
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1 :pvec-pprod 



We will heavily use also the following statement about composition of paravector-field and 
paraproduct operators (see the appendix in [12] for its proof). 

Lemma 3.11. Fix s €]0, 1[. There exist constants C, depending only on s, such that, for all 
ti < and t2 e M, 

\\TxTuv\\^s-i+t-^+t2 < C'(||-^||c= ||^i|lc*i lbllc*2 + 

+ ll^llc«2 \\Txu\\cs~i+H + \\u\\cti \\Txv\\c^-i+t2) , 

and this is still true in the case ti = with \\u\\loo instead of \\u\\(^o. 
Moreover, if s — 1 + ti + t2 > 0, then we have also 

\\Tx R{u,v)\\cs-i+t-,+t2 < C{\\X\\cs \\u\\cti \\v\\ct2 + 

+ ll^llc*2 \\Txu\\cs-i+H + \\u\\cH \\Txv\\cs-i+t2) . 

3.3 Transport and elliptic equations 

System ([1]) is basically a coupling of transport equations of the type 

r dtf + vVf = g, 
I f\t=o = fo ■ 



(T) 



So, we often need to use the following result, which enables us to solve (T) in the Holder spaces 
framework. 



p: transport Proposition 3.12. Let a > (a > —1 if divv = 0). 

Let fo € , g G L-'^([0, T]; C^) and v be a time dependent vector field in Ch{[0,T] x R^) such that 

Vv G Li([0,r];L°°) if a<l, 
Vv G L\[0,T];C^-^) if a>l. 



ist :no-loss-l 



Then equation (T) has a unique solution f in the space ^P|^/^^ C([0, T]; )) P| C^([0, T]; C^). 
Moreover, for all t G [0, T] lye have 

ft 



JO 



dr 



eq: elliptic 



\: ellipticity 



1 : laxmilgram 



with V'{t) :-- 



!|V^;(t)||L- ifa<l. 



\Vv{t)\\c<T-i if a > 1. 

If f = v then, for all a > (a > -1 if divv = 0), estimate 1^ holds with V'{t) := || V/(t) ||loo . 
Finally, we shall make an extensive use of energy estimates for the following elliptic equation: 

(12) - div(aVn) = divF in , 
where a = a{x) is a given suitably smooth bounded function satisfying 

(13) := inf a(x) > 0. 

We shall use the following result based on Lax-Milgram's theorem (see the proof in e.g. |13)). 

Lemma 3.13. For all vector field F with coefficients in L^, there exists a tempered distribution 
U, unique up to constant functions, such that VII G L'^ and equation (|12p is satisfied. In addition, 
we have 

a* ||Vn||i2 < \\F\\l2 . 
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4 Propagation of striated regularity 



s : striated 

Now we are ready to tackle the proof of theorem 12.31 We will carry out it in a standard way: 
first of all we will prove a priori estimates for smooth solutions to ([1]). Then, we will construct a 
sequence of regular approximated solutions. Finally, thanks to the the just proved upper bounds, 
we will get convergence of this sequence to a solution of our initial system, with the required 
properties. 

4.1 A priori estimates 

First of all, we will prove a priori estimates for a smooth solution (p, u, VII) to system ([1]). 



ist:rho L^inf 



!t :Drho_L~inf 



est :u_L'"p 



l:Du_L~iiif 



4.1.1 Estimates for the density and the velocity field 

From first equation of it follows that 

p{t,x) = po {ipr^{x)) , 
so, as the flow ipt is a diffeomorphism over at all fixed time, we have that 
(14) < p* < p{t) < p* . 

Applying the operator di to the same equation, using classical estimates for the transport 
equation and Gronwall's lemma, we get 



< ||V/3o||l°° exp J \\Vu\\l^ dr^ . 



(15) l|Vp(t)||L° 

From the equation for the velocity, instead, we get in the same way 



/* 


vn 


Jo 


p 



dr ; 



LP 



\\u{t)\\LP < \\uo\\lp + 

SO, using (|14p and Holder inequalities, the following estimate holds for some 6 g]0, 1[: 
(16) 



\u 



LP 



< 



\uo\\lp H 

p* 



\VU\\l2 \\s/u\\l-JdT. 



Remark 4.1. Let us observe that, as regularity of the pressure goes like that of the velocity field, 
one can try to estimate directly the norm of the pressure term. Unfortunately, we can't solve 
its (elliptic) equation in this space without assuming a smallness condition on the gradient of the 
density. So, we will prove that VH is in n L°°, which is actually stronger than the previous 
property and requires no other hypothesis on the density term. 

From (|15p it's clear that we need an estimate for the norm for the gradient of the velocity. 
As remarked before, we can't expect to get it from the hypothesis Q G the key will be the 
further assumption of more regularity of the vorticity along the directions given by the family Xq . 
Here we quote also a fundamental lemma, whose proof can be found in p] (chapter 7) for the 
2-dimensional case, in |T2] (section 3) and [18J (again section 3) for the general one. It is the key 
point to get the velocity field to be Lipschitz and it turns out to be immediately useful in the 
sequel. 

Lemma 4.2. Fix e g]0, 1[ and an integer m > N — 1 and take a non- degenerate family Y = 
i<A<m vector-fields over M such that also their divergences are in . 
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Then, for all indices 1 < i,j < N, there exist functions Uij, 6^-^ (with I < k < N, 



1 < X < m) such that, for all (x,^) S M x M , the following equality holds: 



Moreover, the functions in the previous relation could be chosen such that 



and 



m 



- ^' I{Y) 



IV-III 9Ar-10 
I \\\C' 



Now, we can state the stationary estimate which says that the velocity field u is Lipschitz. 
This can be done as in the classical case, because it's based only on the Biot-Savart law. 



p:Du_L~inf Proposition 4.3. Fix e G ]0, 1[ and ^ € ]1, +oo[ and take a non- degenerate family Y = (^A)i<A<m 
of vector-fields over such that also their divergences are still in . 

Then there exists a constant C , depending only on the space dimension N and on the number 
of vector-fields m, such that, for all skew- symmetric matrices with coefficients in L'^ (^Cy, the 
corresponding (by Q ) divergence-free vector-field u satisfies 



est :Du_L~inf 



:i7) 



|Vn||^oo < C 



1 



8(1-8) 



I^IIlcx) log e + 



4.1.2 Estimates for the vorticity 



As in |14) . using the well-known L'^ estimates for transport equation and taking advantage of 
Gronwall's lemma and Holder inequality in Lebesgue spaces, from ([6|) we obtain, for some 7 G 
]0,1[, 



est:Om_L-q (18) < C exp || VnHioodr^ 



X ||^^o||l9 + 



Moreover, of course an analogue estimate holds also for the L°° norm: 



vn||^, \\vu\\l-Jdr 1 . 



est : Om_L~inf 



(19) 



\n 



< C exp 



\Vu\\LoodT I X 



X I ||rio||L°o + 



1 



(p*) Jo 



|Vn||^oc dr 



r:q 



Remark 4.4. Let us fix the index p pertaining to u and let us call q the real number in [2, -|-oo[ 
such that 1/p -\- 1/q = 1/2. From our hypothesis, it's clear that q < q; therefore, thanks to 
Holder and Young inequalities, we have 

linil - < linil'' IIOII^"'' < IIOIl 

ll^'ll^q ^ II'^'IIl? — II liL'JnL"" ' 



sss : est-p 



4.1.3 Estimates for the pressure term 

Now, let us focus on the pressure term: taking the divergence of the second equation of system 
([1]), we discover that it solves the elliptic equation 



eq:Pi (20) 



V p 



div (u • Vii) . 
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est :Pi_L-2 



From this, remembering om' hypothesis and remark \4A\ estimate ([5]) and lemma [3.13| the control 
of norm immeditely follows: 

(21) \ ||Vn||i2 < C \\u\\lp llf^llL^nL- . 

p* 

Moreover, we have that VII belongs also to L°°, and so VII € for all a £ [2, +00]. In 
fact, we are going to show a stronger claim, that is to say VII € C}. Cutting in low and high 
frequencies, we have that 

||vn||ci < ||A_iVn||ci + ||(id -A_i)vn||ci < c(||vn||i2 + ||An||co) . 

Now, from (|20p we obtain 
eq:Lapl-Pi | (22) - AH = V (log p) ■ VU + p div {u ■ Vu) . 

From this relation, the fact that div (n • Vn) = Vn : Vn and the immersion L°° ^ C^, we get 

||An||co < !|An||Loc < ||v(iog/j) • vnii^oo + Hp div (li • Vn)||^oo 

< c(||V/9||^oo livniiioo + p* llVnll^o.) . 

Now, by interpolation we have, for some (3 G ]0, 1[ depending only on the dimension A^, 

livniiioo < c||vn||^_^/2 ||vn||j7'' < cHvniiJ^ ||vn||j7^. 

Thanks to Young's inequality, from this relation and (|2ip one finally gets 
(23) l|Vn||ci < C ((1 + llVpllioo) II^IIlp llS^llL^nL- + P* ||Vu||^„ 



est:Pi_C~l_* 



for some 6 depending only on /3, and so finally on the space dimension N. So our claim is now 
proved. 

Finally, we want to find a bound on the second derivatives of the pressure term. For doing 
this, we will need the striated norm of VII. In fact, passing in Fourier variables and using lemma 
14. 2| for all 1 < i , j < we can write 

Ci^jUiO = a.j(x)iepn(e) + ^6^/(x)(XA(x).e)^fc%). 

k,X 

Applying the inverse Fourier transform and passing to L°° norms, we get 

llv^njl^^ < c(||An||^^ + iiaxvnii^^) . 

Proposition 2.104 of [2] tells us that 

Iia.vn||,„ < f liawnil,, i„g + Mk) . 

Using Bony's paraproduct decomposition to handle the norm in and noticing that the function 
C I—)- ^log(e + k/Q is nondecreasing, we finally get 



est :D-2-Pi 



(24) ||v'n||^^ < c^^iivpii^oo livniici + p* \\Vu\\U + 

+ ||X||c^' llVnlUi log e + 



\\dxS7U\\c. 

ixiir-livnii 
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4.2 A priori estimates for striated regularity 



l:div 



est :div-X 



est : I 



After having established the "classical' estimates, let us now focus on the conservation of striated 
regularity. The most important step lies in finding a priori estimates for the derivations along the 
vector-field X. So, let us now state a lemma which explains the relation between the operators 
dx and div ( • X) (see also remark \2.2\ . 



Lemma 4.5. For every vector-field X with components and divergence in , and every function 
f for some rj € ]0, 1], we have 

||div(/X) - 9x/|lc-i"{...} < cJ\X\\c^\\f\\c.. 

Moreover, the previous inequality is still true in the limit case r] = 0, with \\ ■ \\loo instead of \\ ■ W^o. 

Proof. The claim immediately follows from the identity div (f X) — dxf = /divX and from 
Bony's paraproduct decomposition. □ 

4.2.1 The evolution of the family of vector-fields 

First of all, we want to prove that the family of vector-fields X{t) = (^A(i))i<A<m' where each 
X\{t) is defined by ([7|, still remains non-degenerate for all t, and that each Xx{t) still has 
components and divergence in . Throughout this paragraph we will denote by Y{t) a generic 
element of the family X{t). 

Applying the divergence operator to an easy computation shows us that divl" satisfies 

(a* tt ■ V) divy = 0, 

which immediately implies divy(t) € for all t and 

(25) ||divy(t)||ce < C lldivyollc^exp j \\Vu\\L^dT^ . 

Moreover, starting again from ([8|), we get (for the details, see proposition 4.1 of |12) ) 
{dt + u-V)i A Xxj = *Vn - f A Xx 

from which it follows 

(Vxa) {t,x) = (^""a' Xx^ {0,^i\x)) - 1^ *Vn- (^Vxa) (t, Vr'(V'r(x))) dr , 



This relation gives us 



< 



A Xx {t, x) 



+ 



+ / l|V^x(t-T)|l^. 

Jo 



dr . 



and by Gronwall's lemma one gets 



A Xx {t,x) 



> 



("^a'Xo.a) (Vr^(x)) 



; Jq ||V«||£,oo dr 



Therefore the family still remains non-degenerate at every time t: 
(26) I{X{t)) > I{Xo) exp(^-c \\Vu\\l-o dr^ . 
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est:d X-u 



est:X_C~e 



l:f->Df 



Finally, again from the evolution equation ([8|, it's clear that, to prove that Y{t) is of class C^, 
we need a control on the norm in this space of the term dyu- To get this, we use, as very often 
in the sequel, the paravector-field decomposition: 

dyU = TyU + {dy — Ty) U , 

with (by lemma I3.10p 

\\{dy - Ty) u\\ce < Cj\Y\\ce\\\/u\\Lo^. 
Moreover, for all 1 < i < thanks to ([3]) we can write 



Tyu' = -Y,[dk{-^r^Tyjd,Vl,k - [dk{-^r\Ty,dj 



ik 



Obviously, by lemma [3. 101 again we have 



ik 



< \\Tyn\y-i < \\dyn\\ce-i + c\\y\ 



while for the commutator term we use lemma 13.71 which gives us the following control: 



dk{-A)-\Ty,d^ 



ik 



< C\\Y\ 



101 



So, in the end, from the hypothesis of striated regularity for the vorticity we get that also the 
velocity field u is more regular along the fixed directions and 



(27) 



\\dyu\\cs < c (\\dyn\\ce-i + ||y||c. iivuIIlo 



Moreover, applying proposition 13.121 to (|8]) and using (|27p . (|25p and Gronwall's inequality finally 
give us 



(28) ||y(t)||c^ < C exp(c I^WVuho^dT^ (^\\Yo\\c^ + j\-^fo\\^''\\ 



. dr' 



\dy^\\i2e-i dr 



These estimates having being established, from now on for simplicity we will consider the case 
of only one vector-field X{t): the generalization to the case of a finite family is quite obvious, and 
where the difference is substantial, we will suggest references for the details. 

4.2.2 Striated regularity for the density 

Now, we want to investigate propagation of striated regularity for the density. First of all, let us 
state a stationary lemma. 

Lemma 4.6. Let f be a function in C}. 

(i) If dxf £ and V/ € then one has dx^f G C^^^ and the following inequality holds: 

\\dxVf\\c.-i < C(\\dxf\\c. +J\X\\c^{\\f\\ci + WVfUo 



est:f->Df (29) 



(a) Conversely, if dx^f G ^, then dxf € and one has 
est:Df->f| (30) \\dxf\y < C (J\X\\cs {\\f\\ci + ||V/||i^) + \\dxVf\y-i 
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est:d-T_D 



Proof. (i) Using the paravector- field operator (remember definition 13. 9p . we can write: 

dxVf = [dx - Tx)Vf + TxVf. 

From lemma 13.101 we have that the first term of the previous equality is in C^^^ and 

(31) Widx - Tx) V/llc.-i < c]\X\\c. \\Vf\\L^ . 

Now, we have to estimate the paravector-field term: note that 

TxVf = V{Txf) + [Tx,V]f- 

From the hypothesis, it's obvious that V (Txf) £ C^~^. For the last term, remembering 
that V and Tx are operators of order 1, we can use lemma ISTl and get 



est:T D 



(32) 



[Tx,V]/||c.-i < C||X||c.||/lb 



r:lein f->Df 



est : d_X-rho 



est:d_X-Drho 



Putting together (|3ip . (|32p and the control for ||V (Tx/)||(je-i gives us the first part of the 
lemma. 

(ii) For the second part, we write once again dxf = Txf + {dx — Tx) f ■ 

By definition of the space C^, we know that V/ is bounded: so, the latter term can be 
easily controlled in thanks to lemma [3.101 Now let us define the operator ^' such that, 
in Fourier variables, for all vector-fields v we have 



{^V) (0 



1 



So, noting that the paravector term involves only the high frequencies of /, we can write 

Txf = Txi^Vf) = ^TxVf + [TxM^f. 
Now, applying lemmas 13.101 and 13.71 completes the proof. 

□ 

Remark 4.7. Let us note that, if / G L'* (for some a G [l,-|-oo]) is such that V/ G then 
f G Cl (indeed / G C^'^) and (separating low and high frequencies) 



\f\l 



< C 



L« + liv/i 



L°° / 



Both u and p satisfy such an estimate, respectively with a = p and a = +00. 

Thanks to lemma |4?6| we can equally deal with p or Vp: as the equation for p is very simple, 
we choose to work with it. Keeping in mind that [X{t) , 9f + n • V] = 0, we have 

dtidxp) + u-v{dxp) = 0, 

from which (remember also (|30p ) it immediately follows that 



(33) \\dx{t)P{'t)\\c^ < II^oIIcHIpoIIci + lli^XoVpollc^-i ) exp ( c 
Therefore, keeping in mind also (|28p . one gets also 

(34) \\dxit)'^p{t)\\^. 



L°° 



dr 



< C exp llVnllioo drj x 

X (ibollci ll-'^ollc^ + ll-^XoVpollc^-i + 
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4.2.3 Striated regularity for the pressure term 

In this paragraph we want to show that geometric properties propagates also to the pressure term, 
i.e. we want to prove 9xVn E . 

Again, we use the decomposition 9xVn = Tx (VII) + {dx — Tx) VII. 

As usual, lemma [3.101 gives us 



\\{dx - Tx)vu\\c. < c\\x\\c^ ||v2n| 



Now we use estimate (|24p . the fact that log(e + C) ^ e + C^^"^ and Young's inequality to isolate 
the term ||5xVn||(j£. As 2z < 1 + z^, we have 



\X\ 



ivni 



< C{\\x\\c^ livnjici + ||x| 



ivnii 



and finally we can control ||(i9x — Tx) Vn||(^£ by the quantity 



t:d X-T X-Pi 



(35) C [\\p\\w^.^ W^Wc- l|Vn||ci + \\X\\ce \\Vu\\l^ + \\X\\^c- l|Vn||ci j + - \\dxVU\\ce . 

To deal with the paravector term, we keep in mind that VII = V {—A)^^ (gi + 52), where 
we have set 

gi = —V (log p) ■ Vn and 92 = P div (n • Vu) . 

So it's enough to prove that both TxV ( — A)^^ gi and TxV (— A)^"*^ g2 belong to C^. 
Let us consider first the term 

eq:T_X-Op| (36) TxV {-A)-^ g2 = V{-A)-^Txg2 + [t^, V (-A)-^] 52 • 

From lemma 13.71 one immediately gets that 

< c\\X\\ce \\g2\\co < Cp*\\X\\ce ||Vtx||ioo , 



est:g_2-comm (37) 



Tx,V(-A) 



-1 



92 



C' 



while it's obvious that 



V{-A)-^Txg2 < C \\Txg2\y-1 • 



Now we use Bony's paraproduct decomposition and write 

Txg2 = TxTpdiv{u-Vu) + TxT^w(^u.yu)P + T^i? (p, div (n • Vn)) 
From proposition 13.51 and the equality div (n • Vn) = Vn : Vn, it follows that 



est:g_2-ppl (38) ||7x7div {n-Vn)P||ce-i ^ C II-'^IIl"" H^div («-V«)P||c£: - ||^|[cnip|lci li ; 



and the same estimate holds true for the remainder term TxR {p, div {u ■ Vn)). Lemma 13. IH 
instead, provides a control for HTxTpdiv (n • Vn)||^e_i by (up to multiplication by a constant) 

ll-'^llc^ llplici l|Vn||ioo + ||Vn||ioo \\Txp\\c^'i + \\p\\ci llTxdiv (n • Vn) , 

where ||rx/9|lce-i ^ C'll^llc^ IIpIIc^ proposition 13.51 Now the problem is the control of the C^~^ 
norm of T^div (n • Vn). Writing 



Txdiv (n • Vn) = ^2 TxTg^^^ dju' + TxdiR{u^ , d, 



u 



^ ITxkdkTg^^dju' +d,Tx,dkR{u\dju') - TQ^x^dkR{u\d,u'), 
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by use of lemma 13.111 we can easily see that it's bounded by 

\\u\\qi ||Vn||ioo + 11^11^1 ||TxVu||(^£-i + ||Vn||i;,oo |[rxti||c£ ■ 
Hence, keeping in mind lemmas 13.101 and 14. 6[ we discover 

||Txdiv (u • Vu)||cs-i < C Q|X||ce \\u\\'^i + ||9xti||c^ ll'"llci 
and therefore 

est:g_2-pp2] (39) HTxTpdiv (n • Vn) < C (^||X||c= + llpllci li^llci llt^x-u-Hc-' 

Putting inequahties (|37p . (|38p and (|39p all together, we finally get 

TxV {-Ay^ g2 < C (\\X\\c^ WpWci \\u\\li + \\p\\ci \\u\\ci \\dxu\\cA , 



est:g_2 (40) 



l:d X-F 



for some constant C which depends also on p* and p*. 
Before going on, let us state a simple lemma. 

Lemma 4.8. Fix a e g]0, 1[ and an open interval / C M. 

Let X he a vector-field with divergence in and F : I —^M be a smooth function. 

Then, for all compact set J C / and all p E 1^^'°° valued in J and such that dxp ^ , one 
has that dx{F op) £ and dx^{F o p) € C^^^. Moreover, the following estimates hold: 

\\dx{Fop)\\^, < C\\p\\wi,^ \\dxp\y 
||9xV(F o p)||^,_i < C|1p|1vki,°° (||<9x/5||c= + i|X||cHlp|lH^i. 
for a constant C depending only on F and on the fixed subset J. 

Proof. The first inequality is immediate keeping in mind the identity dx{F op) = F'{p) dxp and 
the estimate 

||i^'(p)||(^e < C ||-^"||ioo(J) WpWc^ < C \\F"\\laa(^J'j • 

For the second one, we write: 

dxV{Fop) = dx {F'{p)Vp) = F'{p)dxVp + F"{p) dxpVp. 

Let us observe that the first term is well-defined in C^"^, and using decomposition in paraproducts 
and remainder operators, we have 



\\F'{p)dxVp\\c^., < C \\F'{p)\\^,^^ l|5xVp||c.-i . 
Now, the thesis immediately follows from lemma l46l 

Let us come back to gi: using the same trick as in (|36p . it's enough to estimate 



□ 



11^x5 



and 



rx,v(-A) 



-1 
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Again, the control of the commutator term follows from lemma 13.71 



est:g_l-comm (41) 



Tx,V(-A) 



-1 



91 



< C i|X||c. Il^illco < — ||X||c. ilVpllioo ||vn||ci . 

' p* 



For the other term, we use again Bony's paraproduct decomposition: 

Tx9i = T>rv(iogp)Vn + rxrvnV(iogp) + rx^(v(iogp),vn) . 
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Thanks to proposition 13.51 we immediately find 

est:g_i-ppi| (42) ||rxTv(iogp) vn||^,_, < c]\x\\c^ ||Vp||l- ||vn|[ci , 



and the same control holds true also for the remainder. Moreover, a direct application of lemma 
13.111 implies 



est:g_i-pp2| (43) ||rxrv(iogp) vn||^,_, < c(||x||cniVp||L- livniici + ||vn|ici ||rxV(iogp)||c.- 

Now, from lemmas l3. 101 and 14.81 we easily get 

||rxV(iogp)||c.-i < c(\\dxp\\c^ \\p\\ + II^II 



1 ' 



Putting this last relation into (|43p and keeping in mind inequalities (|4ip and (|42p . we find 



est:g_i| (44) T^v (-A) ^ 9i < (II^xpHc-' llvn||ci + ||X||c. IIpII^i.oo ||vn||cij , 



where, as before, C may depend also on p* and p^,. 

Therefore, putting p5]) . (00]) and (j44)) together, we finally get 



est:d X-Pi 



(45) ||5xVn||c. < C7(^||p||vi.i,oo ||9xp||c. l|Vn||ci + ||Vn||ci ||X||cHI/o||^i,oo + 

+ ll-^llc^l|vn||ci + ||p|lvi/i,oo||X||c£||u|[^i + IIpIIcj Ihllci ll^^xwllce^ • 

4.2.4 Conservation of striated regularity for the vorticity 

Let us now establish a control on ||9xf^||c'-i- Applying the operator dx to ([6]), we obtain the 
evolution equation for dx^'- 



eq:vort_str| (46) dt (dx^) + n • V {dx^) = dx Vp A VUj - dx ■ Vn) - dx (*Vu • n) . 



The second and third terms of the right-hand side of (|46p can be treated using once again the 
decomposition 

dx {n-Vu + ^Vu ■ Q) = {dx - Tx) [VL-Vu + *Vn -Q) + Tx{^-Vu + *Vn • Q) . 
Lemma l3.10l savs that the operator dx — Tx maps in continuously; as L°° ^ C^, one has 

\\{dx - Tx) {n-Vu + *Vn-Jl)||^,_i < c|'x||c. ||Vu||l- • 

To handle the paravector term, we proceed in the following way. First of all, we note that, as 
divu = 0, we can write 

{n-Vu + ^Vu ■ n).. = J2 {div}" dku^ - dju^ dku!-^ = J] (dk {u^ diu''^ - dk (u' dju''^'^ . 

k k 

So, we have to estimate the norm of terms of the type TxT^yS/u and TxS/R{u, S/u) . For 
the former we apply directly lemma [3. Ill while for the latter we first use the same trick as in (|36p 
and then lemmas 13.111 and 13.71 

ll^xTvuVnllj^^-i < C (\\X\\c^ l|Vn|||oo + ||Vu||loo \\TxVu\\c6-i^ 

\\TxV R{u,'S/u)\\Qe-i < C [WXWc'WuWciWVuWl^^ + ||TxVn||ce-i + ||Vn||Loo HrxuHc^ ] . 
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Hence, from lemmas 13.101 and 14.61 it easily follows 



est : Om-Du 



(47) \\dx {n-S/u + *Vn • Q) \\^^_, < C \\X\\ce \\u\\lp,oo + \\dxu\\cs \\ 

Now, let us analyse the first term of (|46p . It can be written as the sum of three items: 



dx ( ^VpAVU 



^ (dxp) (Vp A vn) + ^ (dxVp) A vn + ^ Vp A (a^vn) . 

p6 pZ pZ 



So, let us consider each one separately and prove that it belongs to the space ^. 
First and third terms are in fact in L°° ^ C^~^, satisfy 



-^{dxp) (VpAvn) 



VpA (axvn) 



< 



< 



c 



(p* 



c 



(p*) 



\\dxp\\c^ IIVpIIl- livniici 
livpii^.. \\dxvu\y . 



Now, let us find a C^~^ control for the second term. Note that it is well-defined, due to the 
fact that both p and VII are in Cl (the product of a function, cr < 0, with a one is not 
even well-defined). With a little abuse of notation (in the end, we have to deal with the sum of 
products of components of the two vector-fields), we write 

{dxVp)VU = T(a^vp)Vn + rvn(9xVp) + R{dxVp,Vn) ; 

remembering proposition 13.51 and the embeddings Cl ^ C^, we get 

||(9xVp) A vn||c.-i < c \\dxVp\\cs-i livniici . 



In the same way, as ||l/p^||^j < ||l/p^||y(/i,oo) we get 



^ {dxVp) A vn 
p^ 



< 



c 



1 + 



ivpll 



L°° 



^ \\dxVp\\ce~i livnii^i 



{p*y V p* 

So, using also lemma l46l we finally obtain, for a constant C depending also on p^ and p*, 



est: wedge (48) 



dx [ ^VpAVU 



< c\^\\x\\c. IIpII^i.oo livniici + llVplli^ \\dxVU\y + 

+ ||p||vyi,oo ||0xV/9||^E-i ||vn|[^i). 



est : d_X-Om 



(49) \\dxm\y 



< 



Therefore, from equation (|46p . classical estimates for transport equation in Holder spaces and 
inequalities (|47p and (|48p . we obtain 

Cexp(^c^*||Vn||Loodr^ (^||axo^^o|lc-i + j\- foW^'-U^^^' x 

X ||X||cH|n||£p,oo -I- ||9x'"||cHl^lkp''^ + 

+ l'^llc- llp||^i,oo llVHlIci + IIVpIIl- ||5xVn||ce + 
-I- ||/?||vi/i,oo ||Vn||(^i \\dxVp\\c^~i] drj . 
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ss : f in-est 



f-est:Pi 



cond-T 1 



f-est:S-u 



4.3 Final estimates 

First of all, thanks to Young's inequality and estimates (|2ip and (|23p . for all 7] € [0, 1] we have 
(50) llVnll^^ llVnll^-" < |jVn||^2nci ^ ((l + l|Vp||ioo) ||J^||i,nL- +P*l|Vn||io 

So, setting 



L{t) := \\uit)\\Lv + \m)\\L'> 



putting (jlSp and (|50p into (|16p . (jlSp and (|19p . for all fixed T > we obtain, in the time interval 
[0,T], an inequality of the form 

L{t) < C exp(^cj^ llVullioodr^ ^L(O) + ^ ||Vii||^^ / -^^('^)^'^) ' 

with constants C, c depending only on A^, e, and p*. Now, if we define 

(51) T := supjt > ^* (e- L{T')dT' ^ |[Vn(r)||2^) dr < 2L(0)| , 

from previous inequality and Gronwall's lemma and applying a standard bootstrap procedure, we 
manage to estimate the norms of the solution on [0, T] in terms of initial data only: 

L{t) < C L{0) and |[/9(t)||y(/i,oo < C ||po||vi/i'°° • 
From this, keeping in mind (|5Up and (|5ip . we also have 

l|Vn||ic„(i2)nLi(ci) < c(i + \\Vpo\\U)l\o). 



Now, let us focus on estimates about striated regularity. First of all, from (|26p we get that 
the family X(t) remains non-degenerate: I(X{t)) > C I{Xq). 

Now, for notation convenience, let us come back to the case of only one vector-field, which 
we keep to call X, and set S{t) := ||5x(t)f^(i)||^£_i- Let us note that the constants C which will 
occur in our estimates depend on the functional norms of the initial data, but also on the time T. 

From (EHl) and (EH) we find 



\\Xit)\y < C(^\\X0\\C^ + ^ S{T)dT^ 

dx(t)u{t)\\ < C(s{t) +J\X{t)\\c.\\Vu{t)\\ 



while (f55|) and give us 

II^xpIIc. < C and ||9xwVp(t)||^,., < C (^1 + f^S{T)dT^ . 

Before going on, let us notice the following fact, which is a direct consequence of the integral 
condition in (1511): for m = 1, 2 we have 



(52) ^* (^J^ S{T')dT'^ \\Vu{T)\\fo. dT < C S{T)dT. 



We will repeatedly use it in what follows. 
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Now, let us focus on SxVII: for convenience, we want to estimate its -Lj(C'^) norm, starting 
from the bound (|45p and the ones we have just found. 
First of all, we have 

\\dxp\\ce l|Vn||ci dT < C 

Jo 

^ ||Vn||ci||X||c.||p||^i,oodr < ||Vn||^i(ci)||X|Uc»(c.) < + ^ 5(r)dr^ . 

Exactly in the same way, using also Jensen's inequality, we get 

^*||x||3.||vn||ci(ir < c||vn||ii(ci)|'x||ic.(c.) < c(^i + 1^ SHr)dT^ , 

while, keeping in mind the definition of the jCP'^ norm (see remark {4.7^ and inequality (|52p . we 
easily find 

J \\p\\w^'°°\\''J'\\cp-°° \\dxu\\cs dr < C^ + J S{T)dT + J \\Vu\\l^ S{t) dr^ . 

Therefore, in the end we get 

||axVn||^i(c.) <c{l + j^ (1 + ||Vu||Loo)5(r)dT + ^ S\T)d^ . 

Finally, let us handle the term S{t): from (09]), we see that we have to control the Ll norm of 
i, defined by (??). First of all, we have 

j \\X\\ce\\u\\\p,oo dr , j \\u\\cv.°o WdxuW^e dr < C ^1 + J S'(t) dr^ : 

we have just analysed the same items multiplied by ||/>||;;i/i,oo, which we controlled by a constant. 
Moreover, one immediately find 



r ll^plli- Pxvniic. dT < C \\dxVU\\ 
Jo 



while the term ||/o||^i,oo ||Vn||(^i already occurred in considering 9xVn, and so it can be 

absorbed in the previous inequality. Finally, we have 

^ \\p\\w^,oo ||Vn||ci \\dxVp\\ce~i dT <c(^ + S{T)d^ . 

Putting all these inequalitites together, in the end we find the control for S{t) on [0, T]: 

S{t) < C (^S{0) + _y (1 + \\'^u\\L^)S{T)dT + ^ S'\T)dT^ . 

Now, suppose that T was chosen so small that, in addition to (|5ip . for all t E [0,T] one has also 

cond-T_2 1 (53) / 5^(r)dr < 2 5(0). 

Jo 
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Then Gronwall's lemma allows us to get the bound 

\\dxit)m\\c^-i < CSiO) yt€[0,T], 

for a constant C depending only on T, N, p, q, e, and p* and on the norms of initial data in 
the relative functional spaces. 

Let us note that this inequality allows us to recover a uniform bound, on [0,T], for ||Vn||2,oo 
and ||Vn||(ji, which we previously controlled only in Lj. 

Remark 4.9. The lifespan T of the solution is essentially determined by conditions (|5ip and 
(j53p . In section [5] we will establish an explicit lower bound for T in terms of the norms of initial 
data only and we will compare it with the classical result in the case of constant density. 

4.4 Proof of the existence of a solution 

After establishing a priori estimates, we want to give the proof of the existence of a solution for 
system ([T]) under our assumptions. 

We will get it in a classical way: first of all, we will construct a sequence of approximate 
solutions of our problem, for which a priori estimates of the previous section hold uniformly, and 
then we will show the convergence of such a sequence to a solution of ([1]). 

Now, we will work with positive times only, but it goes without saying that the same argument 
holds true also for the backward evolution. 

4.4.1 Construction of a sequence of approximate solutions 

For each n £ N, let us define Uq := SnUQ] obviously Uq € L^, and an easy computation shows 
that it belongs also to the space Bp j. for all cr € K and all r G [l,+oo]. Let us notice that 
Cla-^pr ^b°^ ^° particular we have that Uq G fl j., for some fixed s > 1 and r G [1, +oo] 
such that S^ ,. "-^ C°'^ 

Keeping in mind that [S'„, V] = 0, we have that = Sn^o G L"^ H B^l; in particular, from 
dSI) we get Vu^ G L''. 

Now let us take an even radial function 9 G Cq°(]R^), supported in the unitary ball, such that 
< 9 < I and f^^ 9{x) dx = 1, and set 6'„(x) = 9{nx) for all n G N. We define := 9n * Po'- 
it belongs to B^^ ^ and it satisfies the bounds < p^, < Pq < p* ■ 

Moreover, by properties of localisation operators Sn and of we have: 

• Po^ Pq in and ||Vp^||l°° < c||Vpo||L-=; 

• Uq — 7> in the space and 11^0 ||lp ^ c||iio||Lp; 

• VLq ^ VLq \n L't and WVLqWiq < c\\Qo\\li, ||ilollL°° ^ c||r2o||L°°- 

So, for each n, theorem 3 and remark 4 of |14j give us a unique solution of ([1]) such that: 

(i) G C([0,r"];i?^,,), with < p* < p" < p*; 

(ii) n" GC([0,r"];LPnB^^,,), with G C([0, T"]; n B^^); 

(iii) Vn" G C([0,T"];L2) nLi([0,T"];i3^_,,). 

For such a solution, a priori estimates of the previous section hold at every step n. Moreover, 
remembering previous properties about approximated initial data, we can find a control indepen- 
dent of n G N. So, we can find a positive time T < T" for all n G N, such that in [0, T] the 
approximate solutions are all defined and satisfy uniform bounds. 
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4.4.2 Convergence of the sequence of approximate solutions 

To prove convergence of the obtained sequence, we appeal to a compactness argument. Actually, 
we weren't able to apply the classical method used for the homogeneous case, i.e. proving estimates 
in rough spaces as (a > 0): we couldn't solve the elliptic equation for the pressure term in 
this framework. 

We know that C L-([0, T]; C L-([0, T]; L?') and (Vn")„gi, C 

L°°([0,r];L2) and, thanks to a priori estimates, all these sequences are bounded in the respective 
functional spaces. 

Due to the reflexivity of L?' and and seeing L°° as the dual of L^, up to passing to a 
subsequence, we obtain the existence of functions /9, u and VII such that: 

• ^ p in the space L°°([0, T]; 

• u" ^ u in L°°([0,r];LP) and 

• Vn" ^ Vn in L°°([0,r];L2). 

Nevertheless, we are not able to prove that (p, n, VII) is indeed a solution of system ([1]): passing 
to the limit in nonlinear terms requires strong convergence in (even rough) suitable functional 
spaces. So let us argue in a different way and establish strong convergence properties, which will 
be useful also to prove preservation of striated regularity. 

First of all, let us recall that, by construction, Uq — )• uq in and Jig — > JIq ii^ L'^^ ^i^d {Po)n 
is bounded in W'^'°°. So, for a > big enough (for instance, take a = max {A^/p , A^/g}), we have 
that (/Oo)n' ('"o)n' (^o)n bounded in the space C~°. 

Remark 4.10. It goes without saying that the sequences of Mq and still converge in 
moreover, also pQ in this space. Remember that pQ belongs to the space Cl, which coincides 

(see [8] for the proof) with the Zygmund space, i.e. the set of bounded functions / for which there 
exists a constant Zf such that 

|/(x + y) + f{x-y) - 2f{x)\ < Zf \y\ 

for all X, y € M.^ . So, using the symmetry of 0, we can write 

Po(x) - po{x) = ]-n^ [ 9{ny) {po{x + y) + po{x - y) - 2 po{x)) dy ; 

from this identity we get that Pq — )■ po in and hence also in C~°'. 
Now, let us consider the equation for p": 

^tp" = - ■ Vp*^ . 

From a priori estimates we get that (tt")„ is bounded in L°°([0, T]; C^) and (Vp")^^ is bounded in 
the space L°°{[0,T]; L°°); so, from proposition 13. 5 1 one has that the sequence {dtp"')^ is bounded 
in L°°([0, T]; C~"). Therefore (p"')„ is bounded in C'''^([0, T]; C""^), and in particular uniformly 
equicontinuous in the time variable. 

Now, up to multiply by a (/? G I>(M^) (recall theorem 2.94 of 0) and extract a subse- 
quence, Ascoli-Arzela theorem and Cantor diagonal process ensure us that p" — ?> p in the space 
C([0,T];Cr„,"). 

Exactly in the same way, one can show that (p"')„ is bounded in Cfe([0, T] xM^) and it converges 
to p in this space. 

Finally, remembering that p G L°°([0, T]; iy^'°°) (recall the compactness argument), by in- 
terpolation we have convergence also in L°°([0, T]; C^^^^^) for all rj > 0. 
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We repeat the same argument for the velocity field. For all n, we have 

dtv!" = - n" • Vu" - a" VH" , 

where we have set a" := (p"")"^. Let us notice that, as po, Oo '■= {Pq)~^ satisfy the same 
hypothesis and a", satisfy the same equations, they have also the same properties. 
Keeping this fact in mind, let us consider each term separately. 

• Thanks to what we have just said, (a")„ C Cb([0, T] xR^)nL~([0, T]; C|) is bounded; more- 
over, from a priori estimates, we see that also (Vn™)^ is bounded in the space L°°([0, T]; C|). 
Therefore, it follows that the sequence (o" VII")^ is bounded in L°°([0, T]; C^"). 

• In the same way, as C L°°([0, T]; C^) and (Vu")„ C L°^([0, T]; L°°) are both bounded 
sequences, one has that the sequence (u" • Vu")^ is bounded in L°^([0, T]; C~"). 

Therefore, exactly as done for the density, we get that (u")^ is bounded in C'^'"'^([0, T]; C^"), so 
uniformly equicontinuous in the time variable. This fact implies that — )• n in C([0, T]; Cj~°). 

Finally, thanks to uniform bounds and Fatou's property of Besov spaces, we have that u £ 
L~([0,T];C^i) and, by interpolation, that ^ u in C{[Q,T];C]^J') for all > 0. 

So, thanks to strong convergence properties, if we test the equations on a (/? G C^([0, T]; 5(R^)) 
(here we have set S to be the Schwartz class), we can pass to the limit and get that (p, n, VII) is 
indeed a solution of the Euler system ([1]). 

Before going on with the striated regularity, let us establish continuity properties of the solu- 
tions with respect to the time variable. 
First of all, from 

dtp = -u-Vp, 

as u G C([0, T];L°°) (from the convergence properties just stated) and Vp € L°°([0, T]; L°°), we 
obtain that p G C'^'^([0, T]; L°°), and the same holds for a := p~^. 

Remember that u G L°°([0, T]; L^), Vn and a G L°°([0, T]; L°°). Moreover, as VH G 
L~([0,r];L2) nL°°([0,r];L°°), it belongs also to L~([0, T]; L^). So, from the equation 

dtu = — u ■ Vn — a VII , 

we get that dtU G L'^{{0,T]; LP), therefore u G C{[0,T];LP). 

In the same way, from (|6]) we get that G C{[0,T]; L'^), and hence the same holds true for 

Vn. 

Now, using elliptic equation (|20p and keeping in mind the properties just proved for p and a, 
one can see that VII G C([0,T];L^). So, coming back to the previous equation, we discover that 
also dtU belongs to the same space. 

4.4.3 Final checking about striated regularity 

It remains us to prove that also properties of striated regularity are preserved in passing to the 
limit. For doing this, we will follow the outline of the proof in |10] . 

1. Convergence of the flow 

Let ip"" and ^/J be the flows associated respectively to and u; for all fixed G P(M^), by 
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definition we have: 

\^{x){rit,x)-^Pit,x))\ < 

< 



< 



{u^ir, rir, x)) - uir, V(t, x))) | dr 

\ip{x)iu^-u)iT,rir,x))\ + 

+ |(/j(x)n"(r,^"(T,x)) - (y9(x)n"(r,^(r,x))| dr 



|Vu"| 



\^{x) ir-^) {r,x)\dT + 



+ 



ly^n" — ipu\\^oo dr . 



So, from convergence properties stated in previous part, we liave tliat tp"^ tp \n tlie space 
L°°([0,T];/d + L^J. Moreover, it's easy to see tliat 



w^rmLo^ < cexp 



|Vn"||^^ dr 



wliicli tells us that the sequence is bounded in L°°([0, T]; Id + C°'^). Hence, finally 

we discover that tp"^ ^ tp also in the spaces L°°([0, T]-Id + q-") for all ?? > 0. 

2. Regularity of Oxgip 

First of all, let us notice that, by definition, 

dxoi.)r{t,x) = Xnr{t,x)); 

applying proposition 13. 8| we get 



est:X-o-psi (54) 



t \\L°° 



\x. 



t \\c^ 



which implies that {dxo'P'^)n bounded in the space L°°([0, T]; C^). Now we note that, for 
every fixed (/p G P(M^), we have 



the second term is compactly supported, hence it converges in L°° because of what we have 
already proved. So let us focus on the first one: 

dxo {ip^n - dxo {^^) = div (Xo ifi^P^ - iP)) - ^{^^ - iP) divXo ; 

decomposing both terms in paraproduct and remainder and remembering hypothesis over 
Xq, it's easy to see that 



\dx, i^r) - dxo (^V')iic-i <c\\^r-v'' 



IX, 



Therefore, from what we have just proved, OxqiP"' — > dx^'P in L^{\fi^T]]Ci~^); moreover, 
by Fatou's property, one gets that dxo'P G L°°([0, T]; C^) and it verifies estimate (|54p . So, 
by interpolation, convergence occurs also in L°°[\}i,T]]C^~J^) for all > 0. 

3. Regularity of Xf 

Remembering the definition 

Xt{x) := {dxo{x)'^){t,iJt^{x)) and divXf = (IyvXqo^-'^, 

from proposition 13. Sl it immediately follows that Xt and divX^ both belong to . Moreover, 
the same proposition implies that X" — t- X in the space L°°{\fi,T\;C^~J^) for all > 0, and 
the same holds for the divergence. In particular, we have convergence also in L°°([0, T]; L"^^), 
which finally tells us that Xj remains non-degenerate for all t G [0,T], i.e. I{Xt) > c /(Xq). 
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4. Striated regularity for the density and the vorticity 

Let us first prove that regularity of the density with respect to Xf is preserved in time. To 
simphfy the presentation, we will omit the localisation by 93 G 2?(M^): formally, we should 
repeat the same reasoning applied to prove regularity of dx^ip- So, let us consider 

dx^P^ - dxp = div -X)) - />"div(X"-X) + div ((p" - divX 

and prove the convergence in L°°([0, T]; C^^^). Using Bony's paraproduct decomposition, 
it's not difficult to see that first and third terms can be bounded by ||p"||/,oo — X\\Loa + 
ll/o" ~ while second and last terms can be controlled by ||/o"||loo ||div (X" — 

^)llc'^/2 + \\p" ~ p\\l°° ||divX||(^£:/2, for instance. So, from the convergence properties stated 
for and (X")^, we get that dx^^p^ — >■ dxp in the space L°°([0, T]; C^~^), as claimed. 

Moreover, from a priori bounds and Fatou's property of Besov spaces, we have that dxp G 
L°°([0, T]; C^) and so, by interpolation, convergence occurs also in L°^{[0,T];Ci~J^) for all 
7? > 0. 

Now we consider the vorticity term (again, we omit the multiplication by a 2?(M^) function): 



dx^ = div ((X" - X) f)") - f]" div (X" - X) + 

+ div (X® - divX. 

From the convergence properties of we know that 0" — ?> in L°°{[0, T];C^^) for all 

7? > 0, so for instance also for rj = e/2. From this, using again paraproduct decomposition 
as done before, one can prove that ^ dxft in L°°([0,r];Cj '^^). From a prion 

estimates and Fatou's property of Besov spaces again, this implies dx^ G L°°([0, T]; C^~^), 
and moreover convergence remains true (by interpolation) in spaces L^{[0,T];Ci^^~^) for 



all T] > 0. 

So, all the properties linked to striated regularity are now verified, and this concludes the 
proof of the existence part of theorem 12.31 



4.5 Uniqueness 

Uniqueness in theorem 12.31 is an immediate consequence of the following stability result. 
p:stab Proposition 4.11. Let (^p^ and (p^,u^,Vn^) be solutions of system ([1]) with 

< p* < p\ < /9* . 

Let us suppose that dp := p^ - p^ G C{[Q,T]; L"^) and that 6u := - G C^([0, T]; L^). 
Finally, assume that Vp^, Vtt\ Vu^ and ViP all belong to L^{[0,T]; L°^). 
Then, for all t G [0, T], we have the following estimate: 



\\6p{t)\\^, + \\5u{t)\\^, < Ce^^W (IIMO)IL. + IIM 
where we have defined 

I{t) := / (||Vp2|| +||Vni||^^ + ||Vu2|| +||vn2|| )dr, 
Jo 

Proof. From dt6p + • V6p = — 6u ■ Vp^, we immediately get 

\\6p{t)\\L2 < \\Sp{0)\\l2 + f\\Su\\L2 ||Vp2|| dr. 

Jo 
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s : lifespan 



t:life 



life:T 



Moreover, the equation for 5u reads as follows: 

dtSu + ■ V5u = -Su- Vm^ 



V5U vn2 ^ 

H 1-^ op, 



P 



pi p2 



where we have set 611 = 11^ — U^. So, from standard L'P estimates for transport equations, one 
has: 



\\6u{t)\\L2 < \\5um\L2+C f 

Jo 

Now, let us analyse the equation for V^II: 



L2 pU^W^^ + \\V5Ii\\L2 + II Vtf||^^ \\5p\\l2) dT . 



div T" j = div ^- + ■ V(5?i + 5u ■ Vu^^ 

= div - -p^ 5p + 5u- (Vn^ + Vn^) , 
V P J 

where, to get the second equality, we have used the algebraic identity 

div [v ■ Vw) = div {w ■ Vu) + div (v div w) — div (w div v) . 
So, from lemma [3. 131 we obtain 

\\V6U\\l2 < C(||vn2||^^ \\6p\\l2 + \\6u\\l2{\\Vu^\\^^ + \\Vu^^^)) 
and Gronwall's inequality completes the proof of the proposition. 



□ 



Now, let us prove uniqueness part in theorem 12.31 Let (pi,ui, VIli) and (p^, ti^, VII^) satisfy 
system ([T|) with same initial data (/3o,uo), under hypothesis of theorem 12.31 

As 6u{0) = and u G C{[0,T];LP), Vn G C{[0,T];Li), one easily gets that 6u € C^{[0,T]; L"^). 
Moreover, from this fact, observing that also 5p{0) = 0, the equation for 5p tells us that 6p G 
C ( [0 , T] ; ) . Hence proposition KW can be applied and uniqueness immediately follows. 

5 On the lifespan of the solution 

The aim of this section is to establish, in the most accurate way, an explicit lower bound for the 
lifespan of the solution of system ([T]) in terms of initial data only. 
For notation convenience, let us define 



Ln 



I'^^oIIlp + ||f^o||L9nL° 



and 



|VpoI|l° 



Theorem 5.1. Under the hypothesis of theorem \2.3[ the lifespan T of a solution to system ([T]) 
with initial data {po,Uo) is bounded from below, up to multiplication by a constant (depending only 
on the space dimension N, e, p, q, p.^ and p*), by the quantity 



(55) 



min ■! Lq, ll^i 



X Ln log e + 



llf^oll, 



1 + Lo + m 



o\\c% 



5+3 




1 + lll^( 



+ ll^xVpollc^ 



where 6 > 1 is the exponent which occurs in (|23p . 
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est :U 



Proof. Our starting point is subsection 14.31 with the same notations, moreover we define 

e(t) := L{t) log (e+ f|y) , U{t) := \\Vu{t)\\^^ dr , 

Ait) := WVpim^^ , m :=J\Xit)\\c^ , Kit) := ||ax(t)Vp(t)||^,_, . 

It's only matter of repeating previous computations in a more accurate way. 
Let us notice that, from inequality (|17p . for all time t one has 

(56) U'{t) = ||Vn(t)||L- < Ce(t) : 

we will make a broad use of this fact. 

Now, let us define the time Ti := sup{t > [ U{t) < log 2}. Then, on [0,Ti] we have (from 
(HH and ^) 

A{t) < C Aq and |||o(t)||^yi,oo < C ||po||vi/i'°° > 



lif e-est : L 



and so, keeping in mind (IT6)) . (IT8|) . (IT9|) and (l50|) . we get also 

(57) L{t) < c(^L{0) + (l + A^+^^ J^e\T)dT^ 

||Vn||^2nc,i < c(l + A^+i) Q\t). 
In addition, (|26| imphes I{X{t)) > CI{Xq), while from ([28]) and (|27| it follows 

r(t) < c{tq + j'^S{T)dT^ and \\dxu{t)\\c. < C {S{t) + r(t)G(t)). 



Finally, dM]) and ([Ml) entail 



II^xpIIc. < C((l + ^o)ro + i?o) and ||9xV/>||c.-i < C {{l + Aa)V^ + + j'^ S{T)dT^ . 

From the inequalities we've just established, the control of the striated norm of VII immedi- 
ately follows. 

Let us proceed carefully, as done in subsection 14.31 After some simply (even if rough) manip- 
ulations, we get (up to multiplication by constant terms) 



\\p\\wi,o.\\dxp\\eAWnci < {l + 4-^^j(To + Ro)e\t) 

livniici Ixiic. \\pfwi,o. < (i + 4^') (ro + 1^ s{T)dT^ Q^t) . 

Now, thanks to (a -|- 6)^ < C{a'^ + b^) and Jensen's inequality we infer 

Ix||^.||vn||ci < (ro + 1^ s''{T)dT^ (i + 4)e2(t). 

Finally, the fact that ||ti(i)||/:p,oo < Q{t) implies 

||p||iyi.oo^X||c. |h|||p,oo < (1 + ^o) + ^ S{t) dT^ Q^{t) 

\dxu\y < {i + Ao)e{t)(^s{t) + (To + 1^ s{T)dT^e{t)^ 

< (1 + ^o) (1 + To) {S{t) + 1) Q\t) + (1 + ^o) Q\t) f S{T)dT , 

JQ 
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lif e-est :Pi 



life-def :T 



where, in deriving the last bound, we have used also that Q{t) > 1. 
Let us define 

Mo := (l + 4+^) (1 + r|] + Ro) ■ 
as J S < 1 + J S"^, in the end we get 

(58) ||ax(t)Vn(t)||^, <CMo e\t) (^i + s{t) + s\T)dT^ . 

Now let us focus on the striated norm of the vorticity, estimated in (|49p . Analysing each 
term which occurs in the definition (??) of T, we see that first, second and fourth items can be 
bounded by Vn(t)||^^, and the third one is controlled as in (|58p . up to replace Mq with 

Mo := (1 + Ao)Mo. Finally, we have 

iivniici \\dxVp\\c.-i < (1 + Ao) (i + 4) e\t) (^(i + Ao)To + + s{T)dT^ 

< (l + A^-^^) (1 + To + i?o) Q\t) (^1 + ^* S{T)dT^ . 

So, putting all these inequalitites together, we discover that, in [0,Ti], 

S{t) <c{sq + Mo j'^ e2(r) {l + S{t) + 53(r')dr'^ dr^ , 

and by Gronwall's lemma this finally implies 

S{t) < Ce=/oe^Mrf-(^5o + Mo^*e-^o^®'(^')"^'G2(r) (^1 + 53(r')dr') dr) . 

Let us now define T2 as the supremum of the times t > for which both the conditions 

Mo ^ e^{T)dT < 2Lo and Mo ^ ©^(t) (l + 5^(r')dr'^ dr < 2So 

are fulfilled. Note that, as > 1, also (|56p is verified in [0,T2], so T2 < Ti. Therefore, keeping in 
mind (|57p . in [0,T2] one has 

S{t) < CSo, L{t) < CLo and e{t) < CBo, 

because the function (A, cr) 1-^ Alog(e + cr/A) is increasing both in A and a. 

Let us put these bounds in the integral condition defining T2: we discover that T2 is greater 
than or equal to every time t for which 

Mo eg t < 2 Lo and Mo eg t + Mq eg 5;^ - < 2So. 
Therefore, if we define 

(59) T := J<^^{Lo,So} , 

Mo(l + Lo + So)2 

then both the previous relations are fulfilled, for some suitable constant K. Hence, T < T2, and 
the theorem is now proved. □ 
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Remark 5.2. Let us notice that, in the classical case (constant density), the lifespan of a solution 
was controlled from below by 

T,i := C fllOolli.nL- log {e+ nJl^i,"''^" \\ 

(see also |12)). We have just proved that in our case the lifespan is given by (|55p . instead. The 
two lower bounds are quite similar, even if in our case also the initial density comes into play, 
and there are some additional items, basically due to the more complicate analysis of the pressure 
term. 

Remark 5.3. Note also that, in the two dimensional case, the stretching term in the vorticity 
equation disappears. This fact translates, at the level of a priori estimates, into the absence of 
the first two items in the right-hand side of (??). Nevertheless, as we have seen, the analysis of 
Vn produces terms of this kind: for this reason, in dimension N = 2 we weren't able to improve 
the lower bound ([55]) . 

6 "Holder continuous vortex patches" 

First of all, let us prove conservation of conormal regularity. 

Given a compact hypersurface S C of class C^''^, we can always find, in a canonical way, 
a family X of vector-fields such that the inclusion C holds for all r/ G [e, 1 -|- e]. For 
completeness, let us recall the result (see proposition 5.1 of |12)). which turns out to be important 
in the sequel. 

Proposition 6.1. Let T, be a compact hypersurface of class C^'^. 

Then there exists a non-degenerate family of m = N{N + l)/2 vector-fields X C 7^ such that 
Cl C for allr]e [e,l + e]. 

Hence, thanks to theorem 12.31 we propagate striated regularity with respect to this family. 
Finally, in a classical way, from this fact one can recover conormal properties of the solution, and 
so get the thesis of theorem 12.51 (see e.g. |18j . sections 5 and 6, and |12) . section 5, for the details). 

Actually, in the case of space dimension = 2 , 3 (finally, the only relevant ones from the 
physical point of view) one can improve the statement of theorem 12.51 To avoid traps coming 
from differential geometry, let us clarify our work setting. 

In considering a submanifold S C of dimension k and of class C^'^ (for some e > 0), we 
mean that S is a manifold of dimension k endowed with the differential structure inherited from 
its inclusion in R^, and the transition maps are of class C^'^. 

In particular, for all x € S there is an open ball B C M'^ containing x, and a C^'^ local parametriza- 
tion (f : ^ B n Tj with inverse of class C^'^. This is equivalent to require local equations 
H : B -^R'' of class C^'" such that H\Bns = 0. 

Let us explicitly point out that, when we speak about generic submanifolds, we always mean 
submanifolds without boundary. 

Given a local parametrization ip on U := SPl-B, its differential ip^ : TM.^ TU = TS induces, 
in each point x G M'^, a linear isomorphism between the tangent spaces, (p^:^x '■ T^M.^ — ?> T(^(a.)S. 
Moreover, the dependence of this map on the point x G R'^ is of class C^: in coordinates, 99* is 
given by the Jacobian matrix V99. 

Finally, we say that a function / defined on S is (locally) of class (for a > 0) if the 
composition / o : R*^ ^ R is a-Holder continuous for any local parametrization 93. 

Before stating our claim, some preliminary results are in order. Let us start with a very simple 
lemma. 
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l:D->Hold Lemma 6.2. Let f G L°°(R^) such that its gradient is a-Holder continuous for some a > 0. 

Proof. It's obvious using dyadic characterization of Holder spaces and Bernstein's inequalities. □ 

Now, by analogy, one may ask if this property still holds true for a function defined on 
a submanifold, with Holder continuous tangential derivatives. The answer is yes, with some 
additional hypothesis on the submanifold. 



p: Hold-man 



r: Hold-man 



c : Hold-man 



Proposition 6.3. Let S C be a submanifold of dimension k and of class C^'^, for some e > 0. 
Moreover, let us suppose T, to be compact. 

Let us consider a function / : S — > M, bounded on S, such that dxf G for all vector-fields 

X of class tangent to S. 
Then f eC^'^i^:). 

Proof. Let us fix a coordinate set U := B D (for some open ball B C M^) with its C^'^ local 
parametrization (/? : — > {/, and let us define g := f o ip : ^ M. 
Obviously, g G because / G L~(i;). 

Moreover, for all 1 < i < A; let us set <pt.{di) = Xi : then, Xi is obviously of class C^. Hence 
we have dig{x) = Xi(f){{p{x)), i.e. dig at a point x is the derivation Xi applied to the function 
/, and evaluated at the point ip{x). In our notations, we get dig = (dxif) o (p. 

Therefore, from our hypothesis it follows that \/g G C^, and so, by lemma |6^ g G C^'^(M'^). 

In conclusion, we have proved that / composed with any local parametrization if is of class C^'^ 
on M.^. Therefore / G C^'^(S), and, as S is compact, we can bound its Holder norm globally. □ 

Remark 6.4. Let us note that the operator dx depends linearly on the vector-field X. Hence, 
in the hypothesis of the previous proposition it's enough to assume that one can find, locally on 
S, a family {Xi, . . . ,Xk} of linearly independent vector-fields of class such that dxif G C^(E) 
for all 1 < i < /c. 

Corollary 6.5. Let S C be a compact hypersurface of class C^'^ , and let f G C^(M'^). 

iffec'^^', thenf\seC'''{J:). 

Proof. By proposition 16 . 1 1 and non-degeneracy condition, we can find, locally on S, — 1 linearly 
independent vector-fileds Xi . . . X^v-i, defined on the whole and of class C^, which are tangent 
to E and such that div {f Xi) G C%R^) for all 1 < i < iV - 1. 

Moreover, also the divergence of these vector-fields is e-Holder continuous; therefore, using 
also Bony's paraproduct decomposition, we gather that 

dxj = div ifXi) - f divX, G C%R^) V 1 < z < TV - 1 , 



and hence this regularity is preserved if we restrict dxif only to S. 
So, proposition 16.31 and remark [6.41 both imply that G C^'^(S). 



□ 



Now, let us come back to the situation of theorem 12.51 Moreover, let us suppose that the 
hypersurface Sq is also connected: then it separates the whole space into two connected 
components, the first one bounded and the other one unbounded, and whose boundary is exactly 
TiQ. In dimension 2, this is nothing but the Jordan curve theorem, while in the general case > 3 
it's a consequence of the Alexander duality theorem (see e.g. [20], theorem 3.44). For the sake of 
completeness, we will quote the exact statement and its proof in appendix lAl 

So, let us set Dq to be the bounded domain of whose boundary is ODq = Sq and let us 
define D{t) = ipt[DQ). As the flow ij^t is a diffeomorphism for every fixed time t, we have that 
dD{t) = and also the complementary region is transported by ip: D{tY = ipt{DQ). 

Let us denote by xo the characteristic function of a set O. 
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t:v-patches Theorem 6.6. Under hypothesis of theorem \2.5[ suppose also that the initial data can be decom- 
posed in the following way: 

Po{x) = p}){x)xDo{x) + Po{x)xD^,{x) and floix) = ni{x)xDo{x) + flQ{x) xd^{x) , 

with G C^'^po) and % G C^lZ^o)- 

Then, the previous decomposition still holds for the solution at every time t G [0, T] : 



eq:dec_rho (60) p<yt,x) = p\t,x) XD{t){x) + P%t,x) XD(ty ix) 



eq:dec_vort (61) n{t,x) = n\t, x) XDit){x) + ^""{t, x) XDUri^) ■ 



est :L"l_t-Du 



Moreover, Holder continuity in the interior of the domain D{t) is preserved, uniformly on [0,T]: 
at every time t, we have 

p\t) G C^''{D{t)) and n'{t) G C'{D{t)) 

and regularity on D{t) propagates also for the velocity field and the pressure term: u{t) and Vn(t) 
both belong to C^'%D{t)). 

Proof. First of all, let us recall that, by theorem 12. 5| on [0, T] we have 

(62) [ ||Vu(t)||^oo dt < C. 

Jo 

Thanks to first equation of ([TJ, relation (|60p obviously holds, with 

p^'^{t,x) = pn^PT'i^)) ■ 

So, we immediately get that p^{t) belongs to the space C^'^ {D{t)). Let us observe also that a 
decomposition analogous to (|6Up holds also for a = 1/p, and its components a*'*^ have the same 
properties of the corresponding ones of p. 

Now let us handle the vorticity term. We can always decompose the solution in a component 
localised on D(t) and the other one supported on the complementary set, defining 

n\t,x) := n{t,x)xD{t){x) , il^(t,2;) := n{t,x) XD{t)4x) , 

and therefore obtain relation (|6ip . By virtue of this fact, equation (|6|) restricted on the domain 
D{t) reads as follows: 

dtQ' + u ■ vn' = -{n' -Vu + *Vu • + Va' a vn) , 

which gives us the estimate (keep in mind also (|62p ) 



(Ollc. < c + {\\^" ■ + *vu • 11^, + ||va* A vn||^,) d^ . 

We claim that the first term under the integral can be controlled in . As a matter of facts, by 
([3| we know that the velocity field satisfies the elliptic equation 

N 

in D{t), with the boundary condition (by theorem 12.51 and corollarv 16. 5 P u\dD{t) ^ C^'^{dD{t)). 
So (see theorem 8.33 of [19]) we have that ^t|_D(t) G C^''^{D{t)) and the following inequality holds: 

\\u\\ci,^(D{t)) < C (ll^^(i)llL-(D(i)) + h\dD(t)\\ci,s(9n(t)) + ll^llc-(DW))- 
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Let us note that, as pointed out in |19| , a priori the constant C depends on dD{t) through the 
C^'^ norms of its local parametrizations, so finally on exp ^ |[Vn||Loo(iTy However relation (|62p 
allows us to control it uniformy on [0,T]. Therefore, in D{t) one gets the following inequality: 



Iff • Vti+ *V?i-ff 



\C-{D(t)) 



< c m 



\c^-{D{t)) IP|D{t)llci.-(D(t)) , 



which proves our claim. 

Finally, let us handle the pressure term. From what we have proved, Va* is in C^; so 



|Va* A Vn||^, < C ||Va*||^, 



ivni 



However, we want to prove that an improvement of regularity in the interior of D{t) occurs also 
for Vn. In fact, keeping in mind (|22p . 11 satisfies the elliptic equation 



- An = V(log/>^) -Vli + p'Vu: Vu 

in the bounded domain D{t). Now, from what we have proved, the right-hand side obviously 
belongs to C''{D{t)). Moreover, by theorem 12. 5l and corollarv l6.5| we have Vn|g/)(^) € C^'^{dD{t)): 
in particular, as is compact, Il\Q£){t'^ is continuous. Finally, as D{t) is of class C^'^, it satisfies 
the exterior cone condition (see |15) . page 340). So, theorem 6.13 of applies: from it, we 
gather Ii{t) G C^^^{D{t)). Therefore, Vn(t)|£,(t) G C^'=(D(t)) and its norm is bounded by 



\^^\dD{t)\\ci,ei^QD{t)) + Ir «1lc-(D(t)) ll^niici 



ICi.-(D(t)) 



iVnl 



C-(D{i)) ■ 



Putting all these inequalities together and applying Gronwall's lemma, we finally get a control 
for the norm of fi* in the interior of -D(t), and this completes the proof of the theorem. □ 



app : alg_top 



th:piana 



A Appendix — Complements from Algebraic Topology 

Here we want to prove the following theorem, which we used in section (6] For the technical 
definitions, notions and results, we refer to [20J. 

Theorem A.l. For any dimension N > 2, let T, C be a compact, connected hypersurface 
without boundary. 

Then R^\$] has two connected components (say) B and U, one bounded and the other one 
unbounded, whose boundary is just S. 



The previous result implies in particular that S is orientable (see theorem IA.2p . 
Let us note that, if we already assumed this (redundant) hypothesis in theorem lA.ll then the 
proof would be easier (see e.g. |23)). 

The proof we quote here is actually due to A. Lerario. 

Proof. With standard notations, for a submanifold M C and an abelian group G, we denote 
with 

Hk{M;G) and H^{M;G) 

the /c-th reduced homology and cohomology groups of M. with coefficients in G. 

Let us compactify by adding the point at infinity: in this way, we reconduct ourselves to 
work with the A-dimensional sphere = U {oo}. 

Obviously, M^\S and S^\Y. have the same number of connected components. 

By Alexander duality theorem (see theorem 3.44 of [20]) with coefficients in Z2, we have 



~ H 



N-k-l 



(S;Z2 



V yfc > 0. 
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iso : cohomol 



th: orient 



iso:h tilde 



so : conn_comp 



In particular, this is true for A; = 0: 



Now, as E is compact, connected and without boundary, theorem 3.26 of j20) appHes, and 
gives us 

(independentely whether S is orientable or not). In particular, also i/o('S'^\S; Z2) is isomorphic 
to the same group, and this implies that the homology group (not reduced!) 



(63) 



i/o(5''\5];Z2) ~ i^o(5^'\S;Z2 



has rank equal to 2. But the rank of Hq[M] G) is always the number of the connected components 
of M.. Hence, S'^\S has two connected components, A and B. 
Let us suppose that 00 G ^; then 



AU B 



\E = (^\{oo}) U B. 



Now, as N > 2, U := A\{oo} is still connected. 

Hence, U and B are the two connected components of ]R^\E. 
Moreover, it's easy to see (for instance, by stereographic projection with respect to the point 00) 
that U is unbounded, while B is bounded. 

Finally, obviously dB = dU = T.. □ 

As already pointed out, theorem lA. II entails the following fundamental result. Even if it lies 
outside of the topics of the present paper, we decided to quote it to give a more complete and 
detailed picture of the framework we adopted in section |6j 

Again, the proof is due to A. Lerario. 



(for some N > 2) be a compact, connected hypersurface without 



Theorem A. 2. Let E C 

boundary. 

Then E is orientable. 



Proof. The starting point is relation (|63p in the previous proof. Actually, it holds true for any 
submanifold Ai and any abelian group G: 



(64) 



Ho{M;G) ~ Ho{M;G) G. 



Moreover, it is always true that Hq(A4;G) is isomorphic to the direct product of n copies of G, 
where n is the number of connected components of M: 



(65) 



Ho{M;G) ~ G 



(see |20) for the proof of these facts). 

In the previous proof, we established that the rank of //o('S'^\E; Z2) is 2. Then, by (|65p we 
have that it is still 2 if we consider the homology with coefficients in Z: 

rk{Ho{S^\^;Z)) = 2. 

Therefore, keeping in mind (|64p and the Alexander duality theorem, we gather 



Hois'" \^;Z) 



Z 



Now, by theorem 3.26 of [20j, this last condition is equivalent to the fact that E is orientable. □ 
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